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We show that the dark soliton of the Gross-Pitaevskii equation (GPE) that describes the Bose-Einstein con- 
densate (BEC) density of a system of weakly repulsive bosons, also describes that of a system of strongly 
repulsive hard core bosons at half filling. This connection establishes a relationship between the GPE soliton 
and the magnetic soliton of an easy-plane ferromagnet, where the BEC density relates to the square of the in- 
plane magnetization of the system. This mapping between well known solitons in two distinct physical systems 
provides an intuitive understanding of various characteristics of the solitons. 



I. INTRODUCTION 

The Gross-Pitaevskii equation (GPE) has been central to 
explaining various fundamental aspects of the Bose-Einstein 
condensate (BEC) of a system of weakly interacting bosons 
1 1 1. One of the hallmarks of this system is the dark soliton that 
has been theoretically predicted Q as a unidirectional trav- 
eling wave solution for the condensate density in the case of 
weakly repulsive bosons. This soliton, which dies out when its 
speed approaches the speed of sound, has also been observed 
experimentally [3|. The study of solitary wave propagation 
in BEC remains an active frontier with special emphasis on 
non-GPE dynamics[4|, as well as many-body effects such as 
quantum fluctuations and depletion that cannot be described 
by GPE|H|6l. 

The complex order parameter (i.e., the condensate wave 
function) describing a BEC that satisfies the GPE can be iden- 
tified with the bosonic coherent state average [7] of the boson 
annihilation operator in the continuum version of the Bose- 
Hubbard model. Analogously, in a strongly interacting limit 
described by a system of hard-core bosons (HCB) the appro- 
priate order parameter for this system would be the spin coher- 
ent state average [ 8 ] as the system maps to a pseudospin model 
0. In the continuum description of the dynamics of the ex- 
tended lattice Bose-Hubbard model for hard-core bosons with 
nearest neighbor interactions, this order parameter satisfies an 
evolution equation that is different from the GPE ifTUl . We 
will refer to this equation as the HGPE, where the prefix "H" 
has been used to denote its connection to HCB. 

Unlike the GPE which describes the dynamics of a BEC 
with no depletion, the HGPE evolution encodes both the nor- 
mal and condensate fractions in the bose system. In our recent 
study, we showed that the HGPE supports [10] both a dark 
soliton and an antidark[ 1 1 1 soliton (i.e., a bright soliton on a 
pedestal) for the bosonic density. In other words, the dark soli- 
ton of the BEC of weakly repulsive bosons acquires a partner 
that is bright, when the repulsion becomes extremely strong. 
Apart from being bright, the anti-dark soliton is found to be 
of quite a distinct variety compared to the dark soliton that re- 
sembles the GPE soliton. Away from half-filling, i.e., when 



there is a nonvanishing particle-hole imbalance in the back- 
ground, the dark soliton dies out when its propagation speed 
approaches the speed of sound (like the GPE soliton), while 
its brighter partner is found to persist all the way up to the 
sound velocity [ 10 1. 

In contrast, at half-filling, when the particle-hole imbalance 
is zero, the dark and the antidark solitons in the particle den- 
sity profile of the HGPE become mirror images of each other, 
and they both die out at the sound velocity. Interestingly, in the 
corresponding condensate density, the two solitons become 
degenerate and cease to exist as separate entities. In this spe- 
cial case, the condensate density profile is a dark soliton that 
dies out as its speed approaches the speed of sound IfTUl . 

In this paper, we show that for the half-filled case, the lo- 
calized functional form of the solitary wave solution for the 
density of the hard-core bosonic condensate described by the 
HGPE agrees with that of the dark soliton of the GPE, to a 
very good approximation. In other words, if one focuses on 
the condensate fraction of the hard-core bosonic system, this 
strongly repulsive bose gas supports a solitary wave that has 
almost the same profile as the GPE soliton that describes the 
BEC of a weakly repulsive bose gas. 

Furthermore, since the HCB gas in the spin coherent state 
formulation of the extended Bose-Hubbard model with near- 
est neighbor interactions mimics an anisotropic, easy-plane 
ferromagnetic chain, the well known dark soliton of the GPE 
also relates to the magnetic soliton in the easy-plane ferromag- 
net, where it describes the profile of the square of the in-plane 
magnetization. 

In section II, we begin with the many-body Bose- 
Hubbard Hamiltonian that describes the low-energy behavior 
of bosonic atoms in an optical lattice, and briefly outline some 
steps that lead to the order parameter equations in the weakly 
and strongly interacting limits, respectively. In sections III 
and IV, we first derive the evolution equation for the traveling 
waves in the condensate fraction of the atomic cloud obtained 
from the HGPE, and show that to a very good approximation, 
this equation agrees with the corresponding equation derived 
from the GPE. In section V, we discuss the relationship be- 
tween the density solitons of the BEC and the magnetic soli- 
tons of an easy-plane ferromagnet. 
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II. THE BOSE-HUBBARD MODEL, GPE, AND HGPE 

We begin with the extended lattice Bose-Hubbard model in 
d dimensions, whose Hamiltonian is given by 

H = -'Yytb]b J+a +Vn :j nj +a ]+^2 Unj(rij-l)-(p-2t)nj, 
i,a j 

(1) 

where b^ and bj are the creation and annihilation operators for 
a boson at the lattice site j, rij is the number operator, a labels 
nearest-neighbor (nn) separation , t is the nn hopping param- 
eter, U is the on-site repulsion strength, and p is the chemical 
potential. To soften the effect of strong onsite repulsion, we 
add an attractive nn interaction (V > 0) (although our results 
will be valid for repulsive V as well as for V = 0). Such 
a term may mimic certain characteristics of the long range 
dipole-dipole interaction, that has been considered in several 
recent studies 1121 . The term 2tnj is added so that the terms 
involving t reduce to the kinetic energy expression in the con- 
tinuum version of the many-body bosonic Hamiltonian. 

Conventionally, one defines the order parameter for a many 
boson system to be the thermodynamic expectation value of 
the boson field operator. Invoking the concept of a broken 
gauge symmetry allows this expectation value to be nonzero 
below the BEC transition temperature. It has been argued] 7 1 
that this order parameter may be chosen to be the expectation 
value of the boson annihilation operator in the bosonic coher- 
ent state (also known as the harmonic oscillator coherent state 
or Glauber coherent state) representation of the pure quantum 
state. In this description, it can be shown that quantum fluctu- 
ations are absent. 

The Heisenberg equation of motion for the boson annihi- 
lation operator determined from (QJ, after a bosonic coherent 
state averaging yields the following GPE equation for the con- 
densate order parameter ^(r, t) in the continuum descrip- 
tion. 

- (ft 2 /2m)V 2 * fl + U\^ g \H g - pS! g = ihd t ^! g . (2) 

The GPE with U > provides a very useful characteriza- 
tion of the various properties of the BEC of weakly repulsive 
bosons JT]. Here the condensate density p g = |"3/ g | 2 satisfies 
a continuity equation, and the system shows no depletion. 

The HCB limit (U -> oo) of the Bose-Hubbard Hamilto- 
nian ([T]i has emerged as a useful model to describe various 
characteristics of the BEC of a strongly repulsive bose sys- 
tem. The constraint that two bosons cannot occupy the same 
site can be incorporated in the formulation by using field oper- 
ators that anticommute at the same site but commute at differ- 
ent sites, thus satisfying the same algebra as that of a spin-^ 
system. By identifying bj with the spin flip operator ST, along 

with rij — Sj SJ = | — SJ, the system can be mapped to the 
following quantum XXZ Hamiltonian in a magnetic field: 

H S = - Sj • S j+a - g S|S| +Q ] - J2(9 ~ M) (3) 

3,a 3 



Here g = (t — V)d, where d is the spatial dimensionality. 
As we shall show, g = (t—V) > (see below Eq. Hence 
the HCB system maps to a quantum spin- 1/2 system with a 
Heisenberg exchange interaction t, an easy-plane exchange 
anisotropy g and a transverse magnetic field h z = (g — p) 
along the z-direction. 

This mapping to spins suggests that a natural choice for the 
condensate order parameter ^ s of the HCB system is the av- 
erage of the spin flip operator in the spin coherent state repre- 
sentation ff3l . i.e., 

* s =<S+ > . (4) 

In this representation, it can be shown that [ 10 1 the condensate 
density, p s = |* s | 2 =< S~ X S + > is related to total 
particle number density p —< S~S + > as, 

Ps = P(l - P) (5) 

Setting *i> s = ^fpl exp icf> in the continuum description, the 
condensate order parameter ^ s and the particle density p sat- 
isfy the following equations ifTUl 

ih^ s = -^-(l-2p)V 2 f s -V e y s V 2 p+2gpy s -py s (6) 



P=^ L V-[/7(l-p)V0], (7) 

with the identification ta 2 = ti 2 jm and Va 2 = V e . Equation 
(|6]l (which can also be written as coupled equations for p and 
ip) will be referred to as HGPE. Using the asymptotic value 
p —> p° in Eq. ([6j, the chemical potential p is found to be 

p = 2gp°. 

From Eq. |5]), we note that it is important to distinguish 
between the total bosonic particle density p and the conden- 
sate density p s , where p = p s + p&. Here, pd describes the 
depletion, i.e., the normal component of p. 

Linearizion of HGPE by considering small amplitude fluc- 
tuations in the asymptotic value of ^ s leads to a Bogoliubov- 
like spectrum. For small momenta, the spectrum becomes lin- 
ear, leading to the sound velocity c s for the hard-core system 
as 

c s = (2gp° s /m)K (8) 

where p a s is the asymptotic value of the condensate density. 
Since c s must be real, we obtain the condition 
g=(t-V)> 0. 

In the next section, we will derive an equation of motion 
for the condensate density p s that arises from HGPE, for the 
half-filled case, when the asymptotic particle density is half, 
i.e, when the number of particles equals the number of holes 
in the background. 
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HI. HCB CONDENSATE DENSITY EVOLUTION IN THE 
HALF-FILLED CASE OF HGPE 

By seeking unidirectional traveling wave solutions (along 
the x-direction) of the total density p of the form, 



p{z) = po + f(z), 



(9) 



where z — (x — vt)/a and po denotes the constant asymp- 
totic value of the background density, it has been shown ifTUl 
that the HGPE (Eq. |6]l) leads to a nonlinear differential equa- 
tion for f(z) which can be solved analytically to a very good 
approximation, to yield soliton solutions. 

Here we focus on the half-filled case with po = 1/2. For 
this case, the differential equation for / takes the form IfTUl 



4/ 2 [7 2 - f] 



(10) 



Here, 7 2 = 1 — v 2 with v = v/c s . Further, z — (z and £ = 
A/Vl — A 2 , where the microscopic dimensionless parameter 
A = c s /co is the speed of sound c s in the HCB condensate 
measured in units of the zero-point velocity c = h/ma. 

As discussed in |10|, Eq. (IOi provides the following two 
soliton solutions 

f(z) = ±( 7 /2) sech(2 7 z) = ±( 7 /2) sech(z/r s ), (11) 

where the soliton width F s = (27£)~ 1 . Since 7 must be real 
in Eq. ([TTJ, the dimensionless parameter v must satisfy < 
V < 1. 

The two solutions for / in ( fTT| lead to a doublet of localized 
solitons, namely, a dark soliton and an antidark soliton, for the 
density of the bosonic cloud p(z) = (1/2) + f(z). Using this 
in the condensate density expression given in Eq. Q, we get 



(12) 



Writing p s {z) = (1/4) + f s (z), and comparing it with Eq. 



(12 1 shows that the variation f s around the asymptotic con- 
densate density and the variation / around the total density p 
are related as follows: 



f s {-z)=p s {-z)-- A = -f{z) 



(13) 



Differentiating both sides of Eq. ( 13 1 with respect to z and 
using Eq. ( fTO) in the resulting equation, we obtain the follow- 
ing differential equations for the normalized condensate den- 
sity p s — p s /(l/4) and the corresponding normalized varia- 
tion f s = f s / (1/4), respectively. 



(p) 2 = (l-p s ) 2 (p s -v 2 ) 
az 



{ d -^) 2 = m+i 2 ) 

az 



(14) 



(15) 



Unlike the doublet soliton solutions of / given in Eq. (Hi, 
f s has the unique solution 



fs 



-7 2 sech 2 2 7 z = - j 2 sech 2 z/r s 



(16) 




FIG. 1: (color online) Condensate density profile p s (Eq. ( |12} ) of 
HGPE at half filling, showing the behavior of the dark soliton as its 
propagation speed v = v/c 3 increases from to 1. Note that the 
soliton flattens out at v = 1. 

Figure 1 shows the dark soliton profiles for the condensate 
density for various values of v. As the plot shows, the soliton 
which is dark when 5 = 0, dies out as v 1, i.e., as its 
propagation speed approaches the speed of sound. 



IV. CONDENSATE DENSITY EVOLUTION IN THE GPE 

In this section,we show that Eq. (14) , which describes the 
evolution of the condensate density of the HGPE for a half- 
filled system, has the same form as the condensate density 
evolution equation obtained from the GPE. 

Following the earlier studies of solitons in GPEjT), we write 
the condensate wave function ^ g in that system as 



>P° g {A 



iv) 



(17) 



where p° g is the asymptotic value of the background density. 
tp r denotes the real part of Its imaginary part is v = v/c g 
is the soliton propagation speed measured in units of the sound 
speed c g for the GPE. It is given by 



(Up°Jm)*. 



(18) 



In Eq. ( |17) , < v < 1. Further, it can be shown that ip r 
satisfies Hi 
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where 7 2 



V2 



1 



dlp r 

dw 



( 7 2 - = 0, 



(19) 



and to = {x — vt)/£, where the 



healing length £ is given by the well known expression £ 

h/y/2mUp° = h/V2m c g . 



In analogy with the HGPE (see below Eq. ( 10 1), we define 
a dimensionless parameter A = c g /co, which is the sound ve- 
locity in the GPE measured in units of the zero-point velocity. 
This enables us to write £ in terms of A as 



(20) 



£ = a/V2A 



Hence w — v 2 A(x — vt) /a = v2 Az. Using this in Eq. ( 19 1 
and squaring both sides of this equation, we get 



(^) 2 = (i-V, 2 -^ 2 ) 2 



(21) 



where z = Az. 

Now, from Eq. (17i, the normalized condensate density 
p g = Pg/p g for the GPE is 



Pg = + v 2 



(22) 



Differentiating Eq. (22 1 with respect to z, and using Eq. (21 
it is easy to show that p g satisfies the following equation. 



(|f) 2 = (i-p- 9 m-- 2 ). 



(23) 



The above equation describes the spatial evolution of the 
condensate density whose order parameter obeys the GPE, in 
a frame moving with velocity v. This equation is identical 
in form to Eq. ( 14 1 for the half-filled HGPE, and hence the 



condensate density p g for the GPE (which is equal to to the 
particle density) has the same functional form as the conden- 
sate density p s for the HGPE. We would like to emphasize 
that this mapping between the GPE and the HGPE is valid 
provided one considers only the condensate fraction of the 
HCB. Since the half-filling case for particle density in HGPE 
corresponds to quarter-filling for the condensate density, the 
corresponding GPE density for this mapping to be applicable 
is P° g = 1/4. 

The mapping between the GPE and the HGPE discussed 
above provides the relationship that the width of the soliton for 
the condensate density in the GPE and the HGPE are given, 
respectively, by T g = (2-/A)- 1 and T s = (2 7 A/Vl - A 2 )" 1 

ma. 

We would like to note that the relationship between the GPE 
and the HGPE applies only to the condensate density pro- 
files, as the solitons for the corresponding condensate wave 
functions in these two limiting cases will differ in their phase 
jumps. 



V. THE BEC SOLITON AND THE MAGNETIC SOLITON IN 
THE EASY-PLANE FERROMAGNET 



An interesting consequence of the mapping between the 
GPE and the HGPE discussed in the previous section is the 
emergence of a relationship between two distinct physical sys- 
tems, the GPE soliton in a BEC and a magnetic soliton that 
arises from the quantum easy-plane ferromagnet As we 
discuss below, this provides an interesting way to understand 
why there are in general two solitary waves in the bosonic 
density of the HCB system, and why in the special case of 
half-filling, the condensate part of the HCB has only one soli- 
tary wave and why it is dark. 

In magnetism, the study of strict solitons and solitary waves 
in various types of quasi-one-dimensional classical spin sys- 
tems has been a very active field lfl5l for over three decades. 
Strict solitons arise as solutions of completely integrable sys- 
tems which possess a Lax pair. In view of the relatively few 
examples of systems exhibiting solitons and solitary waves, 
any relationship that is found between different physical sys- 
tems is fascinating, since it unveils the universal mathematical 
aspects underlying these nonlinear systems that describe com- 
pletely different physical phenomena. 

One well known example is the gauge equivalence! 16 1 be- 
tween the completely integrable Landau-Lifshitz equation for 
the dynamics of a spin vector in a classical isotropic Heisen- 
berg ferromagnetic chain and the nonlinear Schrodinger equa- 
tion (NLS), which is just the GPE equation (Eq. Q) in one 
dimension, but with the interaction strength U < (i.e., at- 
tractive interaction). The NLS is well known to support bright 
solitons for the density. As another example, the gauge equiv- 
alence between the classical continuum Heisenberg ferromag- 
netic chain with easy plane (easy axis) single-site anisotropy 
and the GPE (Eq. Q) with U > (U < 0) has been shown 
ifTTll . using their respective Lax pairs. The concept of gauge 
equivalence yields certain useful relationships between 
and the spin field S. Typically, p g = \^g\ 2 gets related to 
the x-derivatives of S and S z ifTTl . Further, magnetic soliton 
solutions have been obtained in the continuum description of 
a classical ferromagnetic chain with a single-site anisotropy 
given by gS 2 z lfl8l 

In the present work, the relationship between the Bose Hub- 
bard model and a quantum ferromagnetic Heisenberg spin 
Hamiltonian (|3]l with an exchange anisotropy and a magnetic 
field will be exploited to relate the BEC soliton to the mag- 
netic soliton. Note that the anisotropy is of the easy-plane 
type since g > 0. 

By writing the complex parameter r appearing in the spin 
coherent representation lfl3l as r = tan(0/2) exp(i<f)), with 
< 6 < 7r and < <j> < 2ir, the coherent state average 
< S + > in the quantum spin- 1/2 system can be calculated in 
terms of 9 and <f>. By identifying p s = | < S + > | 2 we obtain 
the following relation between the condensate density p s and 
the classical spin variables: S x , S y , S z , where 9 and <fr appear 
as the polar and azimuthal angles of a classical spin vector S. 
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dark and antidark solitons are therefore mirror images of each 



Ps = {Si + Si) /A = [(1/4) - sly A 



(24) 



In addition, the particle density p, the chemical potential p 
and the particle-hole imbalance variable (1 — 2p°) in the BEC 
system described by HGPE are related to the spin variables of 
the magnetic system as follows. 



p = 1/2- 
p = 2gp° 
g(l-2p°) = h z , 



S z 



.9(1-25°) 



(25) 
(26) 
(27) 



where h z = (g — p) is the magnetic field along the z- 
direction, as seen from the XXZ Hamiltonian ([3]). Of the 
above, the first equation is obtained by taking the expecta- 
tion value of the corresponding hard-core boson operator to 
be spin coherent state average. The last two are obtained by 
using their asymptotic relationships in the expression for p 
(see below Eq. |7]i). 

On taking spin-coherent state averages of the quantum spins 
on the lattice and going to the continuum, the evolution equa- 
tion for the condensate order parameter also describes the evo- 
lution of classical spins, where the asymptotic density controls 
the external transverse field in the spin system. The solutions 
discussed here correspond to the half-filled case where this 
external field is tuned to zero. 



As a consequence of the relationship (25 I, a soliton solu- 
tion of the total particle density p corresponds to a nonlinear 
excitation of the z-component of the spin. On the other hand, 
the soliton of the condensate density describes a nonlinear 
excitation of the in-plane spin, and relates to the square of the 



other, as seen from Eq. ( 1 1 



In contrast, as seen in Eq. (24 1, the condensate density of 
HCB is related to the square of the in plane magnetization, 
(which maps to the GPE density profile with p° = 1/4) result- 
ing in a single soliton both in the GPE and in the condensate 
fraction of the HCB. Such a soliton has to be a dark soliton, 
since the out of plane component of the spin decreases the 
total magnetization of the system. 

In summary, soliton characteristics in the BEC of bosons 
in weakly and strongly interacting regimes are expected to 
be quite distinct. In cold atom laboratories where the inter- 
action between atoms can be varied, one expects significant 
changes in the dynamics as one tunes the scattering lengths of 
the bosonic atoms. Therefore, our result illustrating the simi- 
larities between the solitons in these two different regimes in 
the special half-filled case is important. Furthermore, by ex- 
ploiting the relationship between the HCB and spin systems, 
the connection between the GPE and HGPE solitons paves 
the way for relating them to magnetic solitons in a ferromag- 
netic spin chain with easy-plane exchange anisotropy. These 
studies open a new way to understand and interpret many im- 
portant characteristics of solitons in these systems, including 
the fact that it provides an intuitive understanding of why the 
GPE soliton is dark. 

By designing a quasi-one-dimensional optical lattice which 
incorporates the hard-core boson constraint of no double oc- 
cupancy, and by loading them with bosonic atoms appropri- 
ately so as to simulate a half-filled lattice, it would be of inter- 
est to study soliton propagation and investigate its relationship 
to the GPE soliton. 



in-plane magnetization, due to the relationship ( 24 1. 



As discussed in our earlier paper] 10], the existence of the 
two solitons in HCB has its root in the fact that both parti- 
cles and holes play equal roles in deciding the dynamics of 
the system. It is noteworthy that the XXZ Hamiltonian (|3]l 
provides an alternative means to understand the existence of 
the two solitons. The spin Hamiltonian <|3j contains terms 
with two competing effects, the easy plane anisotropy g that 
tends to align spins in the xy-plane and the the transverse 
field h z that favors a spin alignment along the z-axis. The 
ground state of the system consists of spins aligned on a cone 
that makes an angle 9o with the z-axis where 8q is deter- 
mined by the background density of the atomic cloud given 
by p Q = 1/2 - S° = 1/2 - (1/2) cos O . 

Interestingly, in the half-filled case p° = 1/2, the particle- 
hole imbalance Sq = (1 — 2p°) is zero, and h z vanishes. In 
this case, the ground state consists of spins in the xy-plane. 
Therefore, the solitary wave excitation resulting in a local- 
ized change in the density makes the spins move out of the 
plane. In view of the fact that there is no clear preferred direc- 
tion, disturbances above and below the easy plane are equally 
preferred. In the corresponding particle density in HCB, the 
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